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Abstract. A total of 19 luminous infrared galaxy merg- 
ers, with L m > 2 x 10 11 L Q , for H = 75 kms" 1 Mpc -1 , 
have been observed in the H I line at Nangay and four 
of them were observed at Arecibo as well. Of these 19, 
ten had not been observed before. Six were clearly de- 
tected, one of which for the first time. The objective was 
to statistically sample the Hi gas mass in luminous in- 
frared mergers along a starburst merger sequence where 
the molecular CO gas content is already known. We also 
searched the literature for Hi data and compared these 
with our observations. 

Key words: Galaxies: distances and redshifts - Galaxies: 
general - Galaxies: ISM - Radio lines: galaxies 



1. Introduction 

Most luminous infrared galaxies (LIGs), with infrared lu- 
minosities Lir^2 x 10 11 Lq, for Hq = 75 kms" 1 Mpc -1 , 
are gas- rich, closely interacting/merging systems (e.g., 
van den Broek et al. 1991; Sanders 1992; Murphy et al. 
1996) with high star formation rates and efficiencies. They 
represent the dominant population among objects with 
£boi~2 x 10 11 L Q in the local universe (Sanders & Mirabel 
1996). 

Most of the LIGs' luminosity is radiated by dust in 
the far-infrared. Evidence in favor of a starburst origin 
for the high IR luminosity is clearly mounting (e.g., from 
molecular gas content, Downes & Solomon 1998; Solomon 
et al. 1997; radio continuum imaging, Condon et al. 1991; 
Smith et al. 1998; ISO mid-IR spectroscopy, Genzel et 
al. 1998; and ground based mid-IR and HST/NICMOS 
near-IR imaging, Soifer et al. 2000; Scoville et al. 2000). 
Although an AGN might also contribute to the high IR 
luminosity of a small fraction LIGs (e.g., Sanders et al. 
1988; Veilleux et al. 1997; 1999; Sanders 1999), the ob- 
served abundant supply of dense molecular gas traced by 
high dipole moment molecules like HCN, usually found 
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only in star formation cores (Solomon et al. 1992; Gao 
1996; Gao & Solomon 2000a, b) indicates that almost all 
LIGs are ideal stellar nurseries. 

A clear correlation between the CO (1-0) luminosity 
(a measure of molecular gas mass M(H2)) and the pro- 
jected separation of merger nuclei (an indicator of merg- 
ing stages) has been shown for a sample of 50 LIG mergers 
(Gao & Solomon 1999), indicating that the molecular gas 
content decreases as merging advances. This suggests that 
the molecular gas content is being rapidly depleted due to 
the starbursts as merging progresses. This was the first 
evidence connecting the depletion of molecular gas with 
starbursts in interacting galaxies. The question arises if a 
similar relationship exists between the Hi gas mass and 
the merger interaction stages. 

Only loose correlations were found between the H I gas 
mass and Lm in previous Hi line studies of LIGs sam- 
ples of LIGs (Garwood et al. 1987; Jackson et al. 1989; 
Martin et al. 1991; Mirabel & Sanders 1988; Young et 
al. 1986), unlike the tight correlation found between Lju 
and the CO gas mass. A tighter correlation was found 
(Kennicutt 1998) between the Lpm/area. and M#j/area, 
however. This may imply that the more diffuse H I gas is 
being depleted as well, like the molecular gas component, 
as merging advances. 

The M(H2)/M(Hi) ratio also increases with Lm for 
LIGs (Mirabel & Sanders 1989). However, whether this 
ratio increases or not as the merging progresses through 
different stages has not yet been explored. With new Hi 
data accumulated for a statistically significant sample of 
LIGs arranged along a starburst merger sequence (e.g., 
by the degree of nuclear separation and the merging mor- 
phology compared with numerical simulations for merg- 
ers), we can ultimately examine how the M(H2)/M(Hi) 
ratio changes as the merger-induced starburst proceeds in 
mergers of different stages, since the molecular gas con- 
tents of our sample galaxies are already known. 

Just over a dozen galaxies in the Gao & Solomon 
(1999) CO sample of LIG mergers have 21cm H I line mea- 
surements reported in the literature. In order to enlarge 
the sample available for a statistical study of the relation- 
ship between the Hi gas mass and the evolution of the 
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merging process as the starburst proceeds along a star- 
burst merger sequence, we have conducted new H I obser- 
vations. Here, we present the results of new Nancay and 
Arecibo 21cm observations of 19 LIGs in the sample (for 
the sample selection, see Sect. 2). 

We will explore possible correlations with the global 
Hi properties of LIGs in detail in another paper (Gao 
et al. in preparation). For this, we will use all Hi obser- 
vations available for LIGs in the CO survey of Gao & 
Solomon (1999), including spectra from an unpublished 
Arecibo Hi survey of LIGs (I.F. Mirabel 1997, private 
communication) . 

2. Sample selection luminous infrared galaxies 
along a starburst merger sequence 

Although there have been extensive Hi and CO observa- 
tions of LIGs, including synthesis imaging (e.g., Mirabel & 
Sanders 1988; Martin et al. 1991; Scoville et al. 1991; Gao 
1996; Gao et al. 1997; Lo, Gao & Grucndl 1997; Downes & 
Solomon 1998; Hibbard & Yun 1996, 1999a, 1999b; Gao & 
Solomon 1999; Gao et al. 1999, 2000a, b), there has been 
no systematic attempt to trace the gas properties (both 
atomic and molecular) during the galaxy-galaxy merger 
sequence starting with systems roughly separated by a 
galactic diameter; e.g., just 4 galaxies in the "Toomre Se- 
quence" (Toomre & Toomre 1972; Toomre 1977), repre- 
senting a merger sequence, were observed by Hibbard & 
van Gorkom (1996). 

However, observations of only the atomic gas along 
the "Toomre Sequence" is not sufficient to trace the star- 
burst processes occurring during a galaxy-galaxy merger, 
as stars are forming in molecular clouds which are prob- 
ably embedded in the more diffuse atomic hydrogen gas, 
and H I itself is not a tracer of the star forming gas. More- 
over, the galaxy mergers in the "Toomre Sequence" are 
optically selected and have large variations in their gas 
content and do not easily fit into a sequence of gas con- 
version into stars. 

Although the "Toomre Sequence" nicely reflects the 
dynamical transformation of two spiral galaxies into a 
merged elliptical, no ultraluminous starburst is included in 
this sequence - the four most infrared luminous objects on 
the sequence are NGC 6621/2, NGC 7592, NGC 2623 and 
NGC 3256, with, respectively, log(L /fl /L Q )=11.26, 11.33, 
11.55 and 11.57. Therefore, although the triggering of a 
starburst by the merging of two spiral galaxies is proba- 
bly represented by the "Toomre Sequence" , it poorly sam- 
ples the ultraluminous starburst phase of gas-rich merg- 
ers along the merger sequence. In short, the "Toomre Se- 
quence" is most likely a morphological sequence, and it 
tells us little about the conversion of gas into stars during 
mergers. Indeed, there are many close mergers where the 
progenitors are not gas-rich galaxies and these may never 
become LIGs (e.g., Bushouse et al. 1988; Kennicutt et al. 
1987). 



The ideal sample to trace a merger sequence leading to 
the formation of an ultraluminous starburst would contain 
gas-rich galaxies of various merging stages that initially 
started with roughly comparably, large gas contents. Since 
this is impossible to ascertain, we selected LIGs where the 
2 progenitors can be identified from CCD images. By mea- 
suring the Hi content as a function of the merger sepa- 
ration for a sample with known molecular gas properties, 
the time dependence of the total (Hi+CO) gas proper- 
ties can then be traced statistically in a starburst merger 
sequence. 

50 LIGs with L IR >2x 10 11 L Q , 2"< S sep ^,(£>i + 
D-2)/2 (where S sep is the separation between nuclei and 
D\ and Di are the major axis diameters of the two galax- 
ies) have been selected for the CO study of Gao & Solomon 
(1999), and we already have molecular gas measurements 
for all these galaxies. Here, we aim at establishing a large 
homogeneous sample of LIGs which have Hi measure- 
ments as well. For observation at Nangay, we excluded 
all high redshift (cz > 20,000 kms -1 ) sources, deemed a 
priori too weak in H I for a lOOm-class telescope, and all 
southern sky (<5 < —39°) objects inaccessible from the site, 
thus reducing the sample to 30 LIGs. By further exclud- 
ing sources around the Galactic Center's right ascension, 
where the transit telescope is too heavily oversubscribed 
with high-priority targets, and by avoiding the repetition 
of previously well-detected sources, we ended up with a 
total sample of 19 LIGs from our CO sample (see Table 
1). For all of these we obtained new global Hi spectra. 

3. Observations and data reduction 

3.1. Nangay observations 

The Nangay telescope is a meridian transit-type instru- 
ment of the Kraus/Ohio State design, consisting of a fixed 
spherical mirror, 300 m long and 35 m high, a tiltable flat 
mirror (200x40 m), and a focal carriage moving along a 
90 m long curved rail track, which allows the tracking of a 
source on the celestial equator for about 1 hour. Located 
in the centre of France, it can reach declinations as low as 
-39°. It has an effective collecting area of about 7000 m 2 
(equivalent to a 94-m diameter parabolic dish). Due to 
the elongated mirror geometry, at 21-cm wavelength it 
has a half-power beam width of 316 E-Wx22' N-S for de- 
clinations below 20°, increasing to 34' N-S at (5=65°, the 
highest declination in the present survey (E. Gerard, pri- 
vate comm.; see also Matthews & van Driel 2000). Typical 
system temperatures were ~40 K for our project. 

Our observations were made throughout the period 
December 1998 - October 1999, using a total of about 150 
hours of telescope time. We obtained our observations in 
total power (position-switching) mode using consecutive 
pairs of two- minute on- and two- minute off-source integra- 
tions. Off-source integrations were taken at approximately 
2 m E of the target position. The autocorrelator was di- 
vided into two pairs of cross-polarized receiver banks, each 
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Fig. 1. Mosaic of Digital Sky Survey images of the 19 infrared luminous merger galaxies observed. The dimensions of 
the individual images are 5'x5'. The order of the galaxies is the same as that used in the Tables. 



with 512 channels and a 6.4 MHz bandpass. This yielded 
a channel spacing of 2.64 kms -1 and an effective veloc- 
ity resolution of ~3.3 kms -1 , which was smoothed to 
a channel separation of 13.2 and a velocity resolution of 
15.8 kms -1 during the data reduction, in order to search 
for faint features. The centre frequencies of the two banks 



were set to the known redshifted H I frequency of the tar- 
get; all galaxies had existing optical and CO redshifts. 

We reduced our Hi spectra using the standard DAC 
and SIR spectral line reduction packages available at the 
Nancay site. With this software we subtracted baselines 
(generally third order polynomials) and averaged the two 
receiver polarizations. To convert from units of T S y S to 
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Table 1. Basic optical data for the sample galaxies. 



No. 


Ident 


R.A. 


Dec 


v 

v opt 


err 


B T 


Dxd 


Morph 


Spect 


sep 






(1950.0) 




















hh mm ss.s 


dd mm ss 


km / s 


km /s 




' x ' 






ii 


1 


IC 1623A 


01 05 19.9 


-17 46 29 


6042 


72 


13.70 


1.0 0.8 


Sd- 


Hn 


16.0 




IC 1623B 


01 05 22.0 


-17 46 18 


5840 


60 


15.00 


0.8 0.6 




Hn 




2 


IRAS 02483+4302 


02 48 19.8 


43 02 53 


15475 


105 


16.75 


0.6 0.3 






3.8 


3 


UGC 2369A 


02 51 15.9 


14 46 24 


9463 


152 


15.39 


0.7 0.3 


Sb 




22.5 




UGC 2369B 


02 51 15.9 


14 46 03 


9380 


118 


15.54 


0.7 0.4 








4 


IRAS 03359+1523 


03 35 57.5 


15 23 10 


10626 


45 


16.58 


0.4 0.2 




Hn 


10.0 


5 


IRAS 04232+1436 


04 23 15.2 


14 36 53 


23972 


60 


18.60 


0.3 0.2 






4.6 


6 


IRAS 08572+3915 


08 57 12.9 


39 15 39 


17480 


42 


(14.92) 






L:/Sy2 


5.5 


7 


IRAS 10035+4852 


10 03 35.7 


48 52 23 


19396 


50 


16.13 


0.6 0.4 






8.7 


8 


IC 2545 

ivy £l\J I ' ) 


10 03 53.0 


-33 38 30 


10245 


47 


15.25 


0.65 0.4 






5.0 






10 03 52.0 


-33 38 42 






(16.83) 










9 


TRAS 10565+2448W 

lllJikJ 1UUU <J | jJT:t:U VV 


10 56 35.5 


24 48 43 


(4 3160") 


(811 










8.0 




IRAS 10565+2448NE 


10 56 35.5 


24 48 43 


(12937) 


(37) 


(16 0) 










10 


Mrk 238 


12 59 20.7 


65 16 06 


14816 


171 


15.92 


0.8 0.7 






14.8 


11 


TR AS 1 3001-2339 


13 00 10.5 


-23 39 13 


6417 


130 


14.68 


1.3 0.5 


I0:p 




6.5 


12 


NGC 5256A 


13 36 14.5 


48 31 47 


8316 


55 


(14 1) 


0.5 0.4 


Comp.p 


Sv2 


9.5 




NGC 5256B 


13 36 15.0 


48 31 54 


8419 


48 


13.67 


0.7 0.6 


Comp.p 


L 




13 


Mrk 463 


13 53 39.6 


18 36 57 


(15289) 




(17.75) 






Sy2 


3.9 






13 53 39.8 


18 36 58 


(14990) 




(17.25) 






Syl 




14 


IC 4395 


14 15 06.2 


27 05 18 


10947 


7 


14.84 


1.1 0.8 


S? 


Hn 


6.0 


15 


UGC 9618A 


14 54 47.9 


24 48 28 


10029 


10 


15.30 


0.8 0.7 


Sc 




41.6 




UGC 9618B 


14 54 48.2 


24 49 03 


10157 


53 


15.70 


0.8 0.3 


Sb 


Hn/L 




16 


Mrk 848 


15 16 19.5 


42 55 30 


12049 


131 


15.95 


0.7 0.3 






6.5 


17 


NGC 6090 


16 10 24.1 


52 35 00 


(9062) 


(30) 


(15.0) 


(0.3 0.3) 


Sd p 


Hn 


6.5 






16 10 24.6 


52 35 04 


(8822) 


(30) 


(14.5) 


(0.4 0.2) 


Sd p 


Hn 




18 


IRAS 20010-2352 


20 01 04.2 


-23 52 25 


15194 


69 


17.07 


0.4 0.3 




Sy2 


8.9 


19 


II Zw 96 


20 55 04.6 


16 56 07 


(10630) 


(81) 


(16.41) 






Hn 


11.0 






20 55 05.0 


16 55 58 


(10770) 


(81) 


(15.20) 











Note: the values in brackets are from the NED database, all others are weighted means from the LEDA database. 



flux density in mJy we used the calibration procedure de- 
scribed in Matthews et al. (2001), see also Matthews ct 
al. (1998) and Matthews & van Driel (2000). This pro- 
cedure yields an internal calibration accuracy of about 
±15% near the rest frequency of the 21cm line, and ±25% 
for the highest redshift source in the sample, at V=24,000 
kms -1 (1300 MHz), as less standard calibration observa- 
tions were made at these infrequently observed high ve- 
locities. 

Another consequence of observing rcdshiftcd objects 
down to a frequency of 1300 MHz, far outside the pro- 
tected 1400-1427 MHz frequency band allocated on a pri- 
mary base to the Radio Astronomy Service (e.g., CRAF 
Handbook on Radio Astronomy 1997), is the occurrence of 
strong radio interference (RFI) signals, which were espe- 
cially noticeable in the 10,000-12,000 kms -1 range (see, 
e.g., the remarks on Mrk 848 in Sect. 4.1). 



3.2. Arecibo observations 

During an observing run at Arecibo from 1 to 13 novem- 
ber, 2000, we obtained H I spectra of the 4 sample galaxies 
which could be observed within the declination range of 
the 305-m diameter telescope and the allocated sidereal 
time slots: IRAS 03359+1523, 04232+1436, 10565+2448 
and II Zw 96. Data were taken with the L-narrow re- 
ceiver, set to linear polarisation, with two polarisations, 
4096 channel subcorrelators and nine channel sampling. 
Each subcorrelator had a 12.5 MHz bandpass, resulting 
in a 1.3 kins" 1 resolution. Observations were taken in 3 
minute ON/OFF pairs, followed by a 10 sec ON/OFF cal- 
ibration pair. The data were reduced using the ANALYZ 
software package (Deich 1990). The two polarisations were 
combined, and corrections were applied for the variations 
in the gain and system temperature of the telescope with 
zenith angle. A 21-channel boxcar velocity smoothing was 
applied, resulting in a resolution of 27.0 kms -1 . Baselines 
were then fitted to the data using the GALPAK package 
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4 IRAS 03355+1 =75 Arocibo 



5 IRAS 04232+1436 Arecibo 



9 IRAS 10565 + 7448 Arecibo 



19 IIZW96 Arecibo 



1 °AA 







1. 02x1 of. 04x1 ot.06x1 Of .08x1 ot 10x1 ot.1 2x1 oil4xlO* 2. 25x 1 o' 1 2. 30x 1 0* 2.35x1 A 2.40x1 o" 1 2.45x1 0^ 2.50x1 o^ 1 . 1 5x 1 0^ 1 .20 x 1 o^ ' .25 x 1 0^ 1 .30 x 1 o" 1 1.35x10* 1.40x10* 9.00x 1 ol.50x 1 0?.00x 1 of ,05x 1 of . 1 0x 1 ot 1 5x 1 0*.20x 1 4 



2 IRAS 02483+4302 



4 IRAS 03359+1523 



-20 1: 

5000 5500 



6500 7000 1.45x10" 1.50x10^ 1.55x10^ 1.60x10" 1.65x10" 8.50x10 J 9.00*1 J 9.50x10 J 1.00x10" 



3.50*10 3 1.00x10* '.05x10 4 1.10x10* 1.15x10" 



5 IRAS 04232+1436 




6 IRAS 08572 + 3915 
I I 




7 IRAS 1 0035 + 4852 




8 IC 2545 



2.30x1 o" 2.35x1 o" 2.40x10" 2.45x10^ 2.50x10^ 1.65x10" 1.70x10" 1.75x1CT 1.80x10" 1.85x10" 1.85x10" 1.90x10" 1.95x10" 2.00x10" 2.05x10" 9.50x10° 1.00x10" 1.05x10" 1.10x10" 



9 IRAS 10565 + 2448 



10 Mrk 238 



RAS 13001-233 



.1 ill 





■? NGC 5256 




.20x10" 1.25x10" 1.30x10" 1.35x10* 1.40x10* 1.40x10* 1.45x10* 1.50x10" 1.55x10" 1.60x10* 5500 6000 6500 7000 7500 7500 SOOO 

13 Mrk 463 14 IC 4395 



8500 9000 



» U in \\ 





1.40x10 1.45X10 1.50X10 1.55x10 1.60x10 1.00x10 1.05x10 1.10x10 1.15x10 1.20x10 9.00x10 J 9.50x10 J 1.00x10 1.05x10 1.10x10 1.10*10 1.15x10 '.20x10 1.25x10 1.30x10 



17 KGC 6790 




18 IRAS 20010-2352 



A 



SOOO 8500 9000 9500 1.45x10 1.50x10 1.55x10 1.60x10 1.00x10 1.05x10 1.10x10 1.15x10 1.20x10 

velocity (km/,) velocity (km/,) velocity (km/,) 

Fig. 2. Arecibo (top row) and Nangay 21-cm Hi line spectra. Velocity resolution is 27.0 kms -1 for the 4 Arecibo 
profiles and 15.8 km s _1 for the 19 Nangay spectra, radial heliocentric velocities are according to the optical convention. 
The spectra are shown in the same order as in Fig. 1 and Tabic 1. The line signal in the Nangay Mrk 848 spectrum is 
caused by residual radio interference, see Sect. 4. 



adapted by K. O'Neil, which was also used to correct the 
velocities to the heliocentric system and to determine the 
global H I line parameters. 

4. Results and Discussion 

The 19 Nangay and the 4 Arecibo Hi spectra are shown in 
Fig. 2. Raw, measured radial velocities, Vhi (in the optical 
convention) and their uncertainties, integrated line fluxes, 
I hi, velocity widths at 50% and 20% of peak maximum, 
W50 and W20, and rms noise levels of our spectra were 
determined using the SIR data reduction package for the 
Nangay data and the adapted GALPAK package for the 
Arecibo data (see Matthews et al. 2001, and references 
therein). 

No corrections have been applied to these values for, 
e.g., instrumental resolution or cosmological stretching 
(e.g., Matthews et al. 2001). We estimated the uncertain- 



ties, <Jv HI , in the central H 1 velocities following Fouque et 
al. (1990): 



a VHI = 4i?°- 5 P°> 5 A- 1 [km/ a] 



(1) 



where R is the instrumental resolution in kms -1 , 
Pvk=(W2o-W5o)/2 [km/s] and X is the signal-to-noise ra- 
tio of a spectrum, which we defined as the ratio of the 
peak flux density and the rms noise. According to Fouque 
et al., the uncertainty in the linewidths is 2<tv hi for W50 
and 3<tv hi for W2o- These Hi profile parameters are listed 
in Table 2, together with values found in the literature. 
For non-detections, the estimated upper limits to the in- 
tegrated Hi fluxes are 2a values for flat-topped profiles 
with a linewidth equal to those measured for the CO line 
emission of the objects, either FWHM or half of FWZI 
values (see Table 3, for the CO data references see Gao & 
Solomon 1999). Derived global properties of the 19 sam- 
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Table 2. Basic Hi data - new results and literature values 



No. 


Idcnt 


R.A. Doc 

(1950.0) 
hh mm ss.s dd mm ss 


rms 
T 

mjy 


Ihi 


Vhi 
k m -l 


W 5 o 
kms -1 


W 2 o 
km s -1 


Tel 


Rot 


1 


IC 1623 


01 05 


21.0 


-17 


46 23 


5.51 


4.56 


6053+20 


287 


331 


N 


* 
















3.30 


6078 


244 


380 


N 


Ma91 


2 


IRAS 02483+4302 


02 48 


19.8 


43 


02 53 


4.90 


<2.45 








N 


* 


3 


UGC 2369 


02 51 


15.8 


14 


46 14 


2.43 


<2.19 








N 


* 
















0.48 


9761 


229: 


— 


A 


Ha97 
















abs. 


9406 


— 


274 


A 


MS88 


4 


IRAS 03359+1523 


03 35 


57.9 


15 


23 11 


1.76 


<0.83 








N 


* 














1.11 


abs. 


10725 


189 


— 


A 


* 
















abs. 


10717 


— 


294 


A 


MS88 


5 


IRAS 04232+1436 


04 23 


15.2 


14 


36 53 


3.43 


<2.74 








N 


* 














0.59 


<0.47 








A 


* 


6 


IRAS 08572+3915 


08 57 


12.9 


39 


15 39 


3.11 


<1.68 








N 


* 


7 


IRAS 10035+4852 


10 03 


35.7 


48 


52 23 


4.85 


<2.43 








N 


* 


8 


IC 2545 


10 03 


52.3 


-33 


38 24 


3.25 


2.56 


10257+16 


182 


228 


N 


* 


9 


IRAS 10565+2448 


10 56 


36.2 


24 


48 40 


2.13 


<1.28 








N 


* 














1.17 


abs 


12825: 


— 


— 


A 


















>0.11 


13100 


— 


>137 


A 


MS88 
















abs. 


12700 


— 


646 


A 


MS88 


10 


Mrk 238 


12 59 


20.7 


65 


16 06 


7.74 


<5.42 








N 


* 


11 


IRAS 13001-2339 


13 00 


10.5 


-23 


39 13 


3.43 


<3.09 








N 


* 


12 


NGC 5256 


13 36 


15.0 


48 


31 51 


4.39 


<2.51 








N 


* 


13 


Mrk 463 


13 53 


39.7 


18 


36 57 


5.93 


<3.55 








N 


* 
















0.20 


15230 


— 


200 


A 


Hu87 
















<0.24 








A 


BB83 
















6.0 


14702 






G 


Hc78 


14 


IC 4395 


14 15 


06.2 


27 


05 18 


3.26 


2.54 


10954+18 


308 


340 


N 


* 
















1.30 


10949 


285 


328 


N 


Ma91 
















0.62: 


10946: 


199: 




A 


Fr88 


1 K 


T TO f Qfil S 


14 54 


48.3 


24 


48 42 


1 42 




QQ8D-I- 1 1 


419 






* 
















3.86 


10029 


317 


442 


N 


Ma91 
















>4.66 


9982 




513 


A 


MS88 
















abs. 


9990 




215 


A 


MS88 


16 


Mrk 848 


15 16 


19.5 


42 


55 30 


4.77 










N 




17 


NGC 6090 


16 10 


24.0 


52 


35 04 


3.57 


4.30 


8901+28 


203 


357 


N 


















3.21 


8841 


211 


232 


N 


Ma91 
















3.4 


8871 




348 


G 


HR89 
















5.29 


8873 






G 


Bu87 


18 


IRAS 20010-2352 


20 01 


04.2 


-23 


52 25 


8.35 


<9.19 








N 


* 


19 


II Zw 96 


20 55 


04.8 


16 


56 02 


2.15 


2.26 


10803 + 17 


270 


315 


N 


* 














1.16 


2.55 


10860+9 


189 


412 


A 


















2.34 


10822 


211 


320 


A 


GH93 
















>2.08 


10810 




263 


A 


MS88 
















abs. 


11132 




300 


A 


MS88 



Notes: ':' denotes an uncertain value; values in italics are for absorption line profiles; 

upper limits to Ihi are 2a values for flat-topped profiles with linewidths assumed equal to the CO line values (see Table 3). 
H I references and telescope codes: 



BB83 Bieging & Biermann (1983) 

Fr88 Freudling ct al. (1988) 

Ha97 Haynes ct al. (1997) 

HR89 Huchtmcier & Richtcr (1989) 

Ma91 Martin et al. (1991) 

-k present paper 

A Arecibo N Nangay 



Bu87 Bushouse (1987) 

GH93 Giovanelli & Haynes (1993) 

He78 Heckman et al. (1978) 

Hu87 Hutchings ct al. (1987) 

MS88 Mirabel & Sanders (1988) 



G Green Bank 90m 



pie galaxies are listed in Table 3; the CO line luminosity 
values are from Gao & Solomon (1999). 

Of these 19 luminous infrared mergers, 10 were 
observed here for the first time in the 21cm line 
(IC 2545; IRAS 02483+4302, 04232+1436, 08572+3915, 
10035+4852, 13001-2339 and 20010-2352; Mrk 238 and 
848; NGC 5256 - i.e., Nos. 2, 5, 6, 7, 8, 10, 11, 12, 16 
and 18). For the 9 objects observed previously, our ob- 
servations served to obtain new, independent Hi profile 
parameters, as in many cases there were discrepancies be- 



tween published data (see Table 2 and Sect. 4.1). Of the 6 
objects in which we clearly detected line emission (II Zw 
96; IC 1623, 2545 and 4395; UGC 9618; NGC 6090 - i.e., 
Nos. 1, 8, 14, 15, 17 and 19), one (IC 2545, No. 8) had not 
been detected before. The H I fluxes measured previously 
at Nangay by Martin et al. (1991) for 4 of these 6 ob- 
jects are on average about 50% lower; we prefer to use our 
Nangay flux values, as they are based on recent external 
calibrations (see Sect. 3). 



7 



In all, 8 of the 19 galaxies observed have reliable H I 
emission line detections: the abovementioned 6 we de- 
tected, as well as Mrk 463 and UGC 2369 (i.e., Nos. 
3 and 13). Two others show an Hi line in absorption: 
IRAS 03359+1523 and 10565+2448 (Nos. 4 and 9), and 
for one (Mrk 848, no. 16) no conclusion could be drawn 
due to strong interference. Global Hi profile parameters 
of (merging) galaxies can depend on the telescope used, 
though, as the different beam sizes (ranging from di- 
ameter at Arecibo to 12' for the Green Bank 90m) may 
include extended gas structures, like tails (e.g., Hibbard 
& Yun 1999a), or nearby, confusing galaxies. And even 
observations made with the same instruments may, sur- 
prisingly enough, sometimes show significantly different 
profiles (e.g., for UGC 2369 and II Zw 96). See the notes 
on individual objects in Sect. 4.1. for details. 

The derived Hi masses listed in Table 3 are plotted 
as function of the projected linear separation between 
the two merger components in Figure 3. Concerning the 
adopted total H I masses for galaxies with more than one 
spectrum available: for the 5 galaxies detected by us we 
used our Nangay results, except for II Zw 96, where we 
used the average of all 3 available emission line fluxes; for 
UGC 2369, we used the Arecibo detection by Haynes et 
al. (1997) and for Mrk 463 we used the Arecibo detection 
by Hutchings et al. (1987) - see also Sect. 4.1. 

Six out of the 9 detected objects have a log(M#j) ex- 
ceeding 9.8, while 2 (UGC 2369 and Mrk 463 - i.e., Nos 3 
and 13) have values ~9.3. A similar plot for the H2 masses 
as function of projected separation, showing a clear corre- 
lation, can be found in Gao & Salomon (1999). The rela- 
tion between the H 1 mass and separation will be examined 
further in Gao et al. (in preparation), using all available 
Hi line data, both published and others, for a sample of 
about 30 objects. 

Notes to individual galaxies 

We searched the vicinity of the target objects for nearby 
spiral galaxies which could possibly give rise to confusion 
in those Nangay H 1 profiles where line emission was 
detected. We used the online NED and LEDA databases 



see Webpages http://nedwww.ipac.caltech.edT; and 



http: / /leda.univ-lyonl ,fr / pages- html / single.html , repec- 
tively), in an area of 5^5x30' (a x S) round the pointing 
centre, i.e. about 1.5 times the HPBW, as well as optical 
images extracted from the Digitized Sky Survey (see 
Fig. 1). For quoted NED values the original literature 
reference is given here (in brackets), while the LEDA 
data listed are weighted average values. 

1. IC 1623 Just inside the Nangay search area, but 
outside the area shown in Fig. 1, lies Sb spiral IC 1622 at 
an E-W separation of 2'.5, i.e. 1.4 times the E-W HPBW, 
SW of the center of the IC1623A/B pair; its Vo pt =6343±60 
kms" 1 , i.e. 285 kms" 1 higher than the mean optical ve- 
locity of the pair, it has magnitude -Bt=14.53, diameter 



Z?25=0?7 and an axial ratio of 0.78. Seen the separation in 
the sky and in velocity, we do not expect our Nangay spec- 
trum to be contaminated by this galaxy, unless it has an 
exceptionally large H 1 envelope. Another nearby compan- 
ion galaxy, nearly 5' NE of IC 1623 at only an E-W sepa- 
ration of 1'.5, was detected in the VLA Hi maps (for pre- 
liminary results, see Hibbard & Yun 1996; also for detailed 
maps, check http : / / www. cv . nrao . edu^j hibbard / vv 1 1 4/| ) , 
which can only contribute less than 10% of our Nangay de- 
tected H 1 emission. Our Nangay data showed strong base- 
line curvature in the H polarisation only, and we therefore 
used the unaffected V polarisation data for the data re- 
duction. 

2. IRAS 02483+4302 In the LEDA database, only 
a single object (PGC 90441) is listed, while two are listed 
in the NED database: Anon 0248+43A, a 16.52 mag ob- 
ject at V=15576±46 kms" 1 (Strauss et al. 1992), and 
Anon 0248+43B, a 17.36 mag object at V=15199 kms -1 . 
We used the LEDA entries for PGC 90441 for Table 1. 
Our Nangay data showed stronger baseline curvature in 
the V polarisation, where we therefore fitted a 5 th or- 
der polyniomial to the data. The data were affected by 
strong, narrow radio interference in both H and V polar- 
isation around 15,525 kms -1 ; this strong, spurious sig- 
nal has been blanked out in Figure 1. At full velocity 
resolution, the narrow absorption feature seen at 15,365 

kms -1 seems to be due to weaker interference, as it is 
much too narrow to be an absorption feature originating 
in the merger galaxy. 

3. UGC 2369 Within the Arecibo beam (3f6 HPBW) 
an H 1 absorption line profile (minimum flux density ^-3 
mJy) was reported centered at V=9406 kins" 1 with a 
W 20 width of 274 kms" 1 by Mirabel & Sanders (1988), 
while an emission line (of average flux density ~2 mJy) 
was reported at V=9761 kms" 1 with a W50 width of 
229: kms" 1 by Haynes et al. (1997), who note that the 
measured width of their profile is poor. The mean optical 
velocities of the 2 galaxies are 9463+152 and 9380+118 

kms" 1 (LEDA), respectively, corresponding very well to 
the velocity of the reported absorption line, but not in- 
compatible with the reported emission line velocity either, 
given the uncertainties in the optical velocities. No line 
signal was noted in our Nangay spectrum, with an rms of 
2.4 mJy. Mirabel & Sanders measured a 21 cm continuum 
flux density of 56 mJy. For the total H 1 mass, we used 
the Arecibo detection by Haynes et al. (1997), which may 
be a lower limit, however, given the broad CO lincwidth 
(FWZI~900 kms" 1 , Gao & Solomon (1999)) and possible 
H 1 absorption at lower velocity. 

4. IRAS 03359+1523 In our Nangay spectrum, with 
an rms of 1.7 mJy, no line signal was noted, neither in 
emission nor in absorption. Within the Arecibo beam (3f6 
HPBW) an H 1 absorption line profile (minimum flux den- 
sity ~-5 mJy) was detected by Mirabel & Sanders (1988) 
and in the present survey. The mean velocity of the ab- 
sorption feature, 10721 kms" 1 , is 95 kms -1 higher than 
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Table 3. Global properties of the sample galaxies 



<1l Sep Lir Leo logM^/ assumed linewidth 

Mpc kpc lO n L 10 9 L; M© kms -1 



1 IC 1623 

2 IRAS 02483+4302 

3 UGC 2369 

4 IRAS 03359+1523 

5 IRAS 04232+1436 

6 IRAS 08572+3915 

7 IRAS 10035+4852 

8 IC 2545 

9 IRAS 10565+2448 

10 Mrk 238 

11 IRAS 13001-2339 

12 NGC 5256 

13 Mrk 463 

14 IC 4395 

15 UGC 9618 

16 Mrk 848 

17 NGC 6090 

18 IRAS 20010-2352 

19 II Zw 96 



81.7 


6.1 


4.7 


213.6 


3.6 


6.2 


127.8 


13.1 


3.9 


144.5 


6.5 


3.3 


327.1 


6.3 


11.1 


236.8 


4.6 


11.9 


263.3 


9.8 


9.3 


137.6 


3.1 


4.5 


174.2 


6.2 


9.6 


204.0 


13.2 


2.5 


86.4 


2.6 


2.4 


110.6 


4.8 


3.1 


204.4 


3.5 


5.3 


147.9 


4.0 


2.2 


137.0 


25.8 


4.1 


162.2 


4.8 


7.2 


117.6 


3.5 


3.0 


206.4 


8.1 


4.7 


149.9 


7.4 


6.6 



10.7 


9.86 




2.9 


<10.42 


250 


7.2 


9.27 




6.9 


abs. 


235 


9.0 


<10.08 


400 


1.7 


<10.35 


270 


7.0 


<10.60 


250 


2.9 


10.06 




5.8 


abs. 




6.8 


<10.73 


350 


2.6 


<9.73 


450 


5.7 


<9.86 


400 


2.8 


9.29 




4.0 


10.12 




17.0 


10.38 




7.0 




300 


5.0 


10.15 




6.4 


<10.96 


550 


6.0 


10.09 





Notes: upper limits to logAfc/ are 2a values for flat-topped profiles with assumed 
linewidths equal to the measured CO line values; abs denotes absorption line spectra; 
Li is in units of K km s -1 pc~ 2 



the optical systemic velocity, which has an uncertainty of 
45 kins" 1 . Mirabel & Sanders measured a 21 cm contin- 
uum flux density of 27 mJy. 

5. IRAS 04232+1436 Both our Nancay and Arecibo 
spectra show a quite narrow absorption line (FWHM=90 

kms -1 at Arecibo, centered at V=23567 kms -1 ), reach- 
ing a minimum flux density of about 5 to 6 mJy. These are 
likely due to interference, however, as the mean velocities 
of these features are about 270 kms -1 different, and the 
measurements were made at 1308 MHz, far outside the 
protected 1400-1427 MHz band. The Arecibo spectrum 
also shows a broader (H / 20 =290 kms -1 ) absorption fea- 
ture at the -2 mJy level, too faint for clear detection in 
Nancay, which may also be due to interference. In conclu- 
sion, the data do not provide sufficient evidence for the 
detection of an H I line. 

6. IRAS 08572+3915 In the LEDA database, only a 
single object (PGC 25295) is listed, while two are listed in 
the NED database (IRAS 08572+3915NW and SE), based 
on near-infrared imaging observations (Zhou et al. 1993); 
as positions only are given for these two components in 
NED, we have used the LEDA entries for Table 1. 

7. IRAS 10035+4852 Very close, interacting galaxy 
pair. In the LEDA database, only a single object (PGC 
29385) is listed, while two are listed in the NED database: 
the E and W objects has optical rcdshifts of, respectively, 



19371+42 and 19421+53 kms -1 (Strauss et al. 1992). We 
have used the LEDA entries for Table 1. 

8. IC 2545 Very close, interacting galaxy pair. In 
the LEDA database, only a single object (PGC 29334) is 
listed, while two are listed in the NED database, based 
on Lauberts & Valentijn (1989). Besides positions, for the 
NW object a radial velocity (10267+86 kms -1 ) and a B 
magnitude (15.27) are listed, and for the SE object only 
a B magnitude (16.83). We have used the LEDA entries 
for Table 1. Detected for the first time in Hi. 

9. IRAS 10565+2448 Two distinct, interacting 
galaxies can be seen on the Digital Sky Survey. In the 
LEDA database, only a single object (PGC 33083) is 
listed, while two are listed in the NED database; we used 
the latter for Table 1. The optical redshift references are 
Kim et al. (1995) for the W object and Strauss et al. 
(1992) for the NE object. Within the Arecibo beam (3'.6 
HPBW) an Hi line profile showing both emission and 
much stronger absorption (maximum/minimum flux den- 
sity ~+l/-5 mJy, respectively) was reported by Mirabel 
& Sanders (1988) - see the two entries in Table 2, where 
the emission line values are lower limits - who measured 
a 21 cm continuum flux density of 44 mJy. Our Arecibo 
observations show some structure in the base line, as they 
were taken after sunrise. The deepest and broadest nega- 
tive feature in our data corresponds well to the one noted 
by Mirabel & Sanders. We lack the sensitivity, however, to 
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Fig. 3. Total H I gas mass plotted as function of the projected linear separation between the two components of the 
16 mergers observed which show line emission or for which an upper limit to the line flux could be estimated; objects 
which show an absorption line only were excluded. Filled circles denote detections and triangles upper limits. Upper 
limits to the Hi masses are 2a values for flat-topped profiles with assumed widths equal to the measured CO line 
width of the objects (see Table 3). 



verify their extremely faint (^-0.5 mJy), very broad ab- 
sorption feature extending down to about 12,300 kms -1 . 
We see no sign of their emission feature adjacent (at higher 
velocity) to the deeper absorption line. 

11. IRAS 13001-2339 = ESO 507-G070 Our 

Nangay data, obtained over 12 days, showed baseline cur- 
vature in both H and V polarisation, and we therefore 
fitted a 5 th order polyniomial to the data. 

12. NGC 5256 Both the LEDA and NED databases 
have entries for 2 objects. In the LEDA database, a B 
magnitude is given for one object only (13.67), while in 
NED the same B magnitude (14.1) given for both objects. 

13. Mrk 463 The Arecibo detection of Mrk 463 
(Hutchings et al. 1987) at 0.2 Jy kms -1 is at the 4ct 
level, with a peak flux density of 1.6 mJy; the authors 
note that the profile was detected in both observing runs. 
Hcckman et al. (1978) note that their Green Bank flux 
measurement is relatively uncertain, 6.0±2.3 Jy kins" 1 , 
and that it would imply a very high Mhi /Lb ratio for Mrk 
463 if all gas detected resides in this galaxy. The reported 
Green Bank flux is 30 times higher than the Arecibo value, 
and the central velocities of these detections seem incom- 
patible: 14702 and 15230 kms" 1 , respectively (AV=528 

kms -1 ), for a W20 profile width of 200 kms -1 measured 
at Arecibo. Inside the Green Bank HPBW but well outside 
the Arecibo beam, the DSS shows only one other galaxy 
that might be a source of confusion in the Green Bank 



Hi spectrum: CGCG 103-016, 6'.1 to the SE of the tar- 
get galaxy, a Bt 15.5 object of O'A diameter and unknown 
rcdshift. If it were the sole source of Hi emission within 
the Green Bank beam, it would have an Mhi /Lb ratio of 
1.4 Mq/L Qi b, which is about 4 times the mean ratio for 
an Sc spiral (Roberts & Haynes 1994), but not unheard 
of for late-type galaxies. The morphology of CGCG 103- 
016 is hard to discern, however, as a relatively bright star 
is situated close to its centre. In conclusion, CGCG 103- 
016 may be the cause of the discrepant Green Bank flux. 
For the total H 1 mass, we used the Arecibo detection by 
Hutchings et al. (1987). 

14. IC 4395 Freudling et al. (1988) note that the pa- 
rameters of their Arecibo H 1 profile are uncertain, as it 
was detected at the edge of the band. Considerable flux 
may well have been missed within the Arecibo band, given 
the fact that our Nangay integrated line flux is consider- 
ably (factor 4) larger than the Arecibo flux and that our 

linewidth of 308 kms -1 is 1.5 times larger than the 
Arecibo value. Our Nancay data showed strong baseline 
curvature in the V polarisation, and we therefore used only 
the unaffected H polarisation data for the data reduction. 

15. UGC 9618 All three available spectra (one from 
Arecibo and 2 from Nangay, see Table 2) show a clear 
emission line with a pronounced central depression. Only 
for the Arecibo spectrum, Mirabel & Sanders (1988) inter- 
preted this as an absorption feature, see their two entries 
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in Table 2; they report a 21 cm continuum flux density of 
109 mJy. 

16. Mrk 848 The Nangay data were heavily affected 
by radio interference, particularly in the H but also in the 
V polarisation, within the range where the H I line emis- 
sion was expected. Though we used only the least affected 
data for the spectrum shown in Fig. 2, the remaining emis- 
sion features appear to be due to RFI. We therefore do not 
consider this galaxy as detected in our survey, but we can- 
not give a proper upper limit to its integrated line flux, 
due to the remnant RFI occuring at its redshift. 

17. NGC 6090 In the LEDA database, only a single 
object (PGC 57437) is listed, while two are listed in the 
NED database. We used the NED entries for Table 1; both 
optical redshifts are from Kim et al. (1995). 

19. II Zw 96 In the LEDA database, only a single 
object (PGC 65779) is listed, while four are listed in the 
NED database; we used the NED entries of the first two 
for Table 1. Near- infrared imaging (Goldader et al. 1997) 
revealed 4 sources; one (I s * in Table 1) seems to be a rel- 
atively quiescent spiral galaxy, another (2 nd in Table 1) 
may be a galaxy nucleus, since it shows hints of a bar 
as well as powerful starburst activity, and the other two 
are most likely highly obscured, luminous starburst re- 
gions. Huge molecular gas concentrations were discovered 
recently around these two highly obscured starburst re- 
gions by Gao et al. (2000a). All published Arccibo, Green 
Bank and Nangay spectra show a clear emission feature. 
Only Mirabel & Sanders (1988) reported an absorption 
feature (minimum flux density -2.5 mJy) in their Arecibo 
profile, adjacent to the emission line (see Table 2), which 
does not occur in the Arecibo spectra taken by us or by 
Giovanelli & Haynes (1993); we therefore consider it to be 
spurious. 
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